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I. INTRODUCTION 

The large time behavior of a system quenched at or below the critical point is characterized by scale invariance . 
For the equal time structure factor one has 

C{k,t) = L"{t)F{kL{t)) (1) 

where 

m^ti (2) 

is the characteristic length growing with time according to a power law and F{x) is the scaling function. The physics 
behind (|l|) is quite simple and is basically due to the degeneracy of the low temperature state. After the quench an 
exponentially fast process takes place leading to local equilibrium. If multiple choice is available, correlated regions 
of the possible low temperature phases are formed. From that point onward equilibration proceeds through the 
coarsening of these correlated regions, whose characteristic size L(t) grows according to (^). The difference between 
quenches to Tc and below Tc is that in the first case the correlated regions are fractal (Appendix I), while in the second 
one are compact. Apart from this, in both cases the equilibration process becomes slow (if the system is infinite, 
equilibrium is never reached) and after domains of the ordered phases have formed, scaling behavior occurs since the 
residual time dependence is confined in the typical size L(t) of the correlated regions. 

The whole time evolution can be divided into a preasymptotic and an asymptotic regime, with a smooth transition 
between the two. The asymptotic regime displays universality and it is controlled by a fixed point structure. The 
universality classes are determined by features like the presence or absence of a conservation law, the number N of 
components of the order parameter, the dimensionality d of space and the final temperature Tp of the quench. More 
precisely, on the temperature axis there is an unstable fixed point at the critical temperature Tc and an attractive 
fixed point at Tp = 0. For the exponent a one has 

r2-^, forT^^T, 
" ~ \ d, for Tf < Tc ^ ' 

where rj is the usual exponent of the static critical phenomena. The exponent z coincides with the exponent of the 
dynamical critical phenomena for Tp = Tc- Instead, for any final temperature below Tc, z = 2 for non conserved 
order parameter (NCOP), while for conserved order parameter (COP) z = 3 when = 1 and z = 4 when N > 1. 
The scaling function F{x) also displays universal features and it is sensitive to the space dimensionality through the 
presence (N < d) or absence (TV > d) of localized topological defects. By contrast, in the preasymptotic regime the 
evolution of the system is not universal, as it depends on the initial conditions of the quench and on the actual value 
of the final temperature. 

A complete theory of the process then should derive the scaling behavior from the basic equation of motion for the 
order parameter and should be able to describe how the relatively simple universal asymptotic regime emerges out 
of the complexity of the preasymptotic regime. Ideally, one would like to have a manageable reference theory which 
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accounts, at least qualitatively, for the basic features of the process and a systematic procedure for the computation 
of the corrections [2|. A scheme of this type is available for quenches to Tc, where, despite the difficulty due to the 
lack of time translational invariance, the field theoretical machinery developed for critical phenomena is to a large 
extent applicable Instead, for quenches below Tc the present status of theoretical understanding is far from this 
standard. What we have in this case is the linear theory Q for the very early stage of the process, which applies only 
when initial conditions are so small that it is actually justified to employ a linear approximation, and ad hoc late stage 
theories ||^. Although these late stage theories have had much success in the computation of the scaling functions, 
yet are based on uncontrolled approximations. Furthermore, the late stage theories do not connect to the early stage 
theory, if this is available at all. So, there is no theoretical understandingof the complex phenomenology arising at the 
breakdown of the early stage theory and leading to the onset of scaling |g] . Proposals for the systematic improvement 
of the late stage theory have been put forward |^] but, as of now, a first principles theory of phase-ordering processes 
is out of reach. 

In this theoretical landscape a special position is occupied by the 1/A'^-expansion. As applied to critical phenomena 
this technique provides a very clear instance of what is to be understood for a systematic theory: there is a lowest order 
analytically tractable approximation (the large- A'^ model) which captures the basic physics and there is an expansion 
parameter (1/A^) which allows for the systematic computation of the corrections. The scheme applies succesfully also 
to quenches to the critical point and, at first sight, it would seem to be applicable as well to the phase-ordering 
processes. Indeed, in the large-A^ model one can solve exactly Jp| for the structure factor and one finds that the 
standard scaling form (|l|) is obeyed for large time with NCOP In particular one finds z = 2, and a is given by 
(3) with ?7 = 0. The scaling functions can also be found explicitely It is to be stressed that in the solution of the 
model with N = oo there are no ad hoc hypotheses and the above outlined picture of the asymptotic behavior with 
scaling, universality and temperature fixed points is derived from the solution of the equation of motion. 

However, when the model is solved with COP although the form (|l|) is obeyed with a — 2 and z = 4 for 

Tp = Tc, for the quenches to Tp < Tc the more general multiscaling form 

C{k,t)^[L{kraL)^r(-^F{x) (4) 
is found, where L{t) ~ t^^^, km{t) is the peak wave vector and x = k/km- The exponent a{x) is given by 

a{x) = g + Qtpix) (5) 

with 



2, for Q<Tp <Tc 
0, for Tp = Q 



(6) 



d, for Tf = 
Furthermore, when < Tp < Tc the function (p{x) in is given by 

?/'(a;), for x < x* 



J d - 2, for < Tj. < Te , . 



with X* — \/2 and 

ij{x) = 1 - {1 - x'^f (9) 
while (p{x) = ipix) for all values of x when Tp = 0. Finally, 



Leaving aside for the moment the apparent formal complication of |1C|) , the important feature which is immedi- 
ately evident is that, contrary to (||), now there are three distinct asymptotic behaviors for Tp = Tc, for < Tp < Tc 
and for Tp ~ 0. For Tp = Tc the structure factor obeys standard scaling with a = 2 as in the NCOP case. Instead 
for Tp < Tc the exponent a depends on x (Fig. |^) and the scaling form (^ involves two lengths, k^[t) and L{t), 
which differ by a logarithmic factor |^ 
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(fc„L)4 = log L'^ + (2 - d) log(fc„L). 



(11) 



The functional form of a{x) is different for < Tp < Tc and for Tp — 0. This means that Tp — Tc and Tp — are 
both unstable fixed points and in between there is a new line of fixed points for < Tp < Tc- The temperature below 
the critical point is no more an irrelevant variable. 

If the 1/N expansion were a good systematic theory, the 1/N corrections ought to produce only minor quantitative 
changes on the picture outlined above. However, this expectation has not been fulfilled by the work of Bray and 
Humayun (BH) who found that for quenches to Tp — and COP standard scaling of the form (|l|) is restored 
in systems with any finite value of N. The same result is very likely to apply also to the quenches to < Tp < Tc- 
Therefore, the main qualitative features emerging in lowest order, like the multiscaling behavior and the relevance of 
the temperature fluctuations, are expected to be a peculiarity of the case with N strictly infinite, disappearing as soon 
as higher order corrections are taken into account. In other words, the limits N ^ oo and i — *■ c» do not commute 
for quenches to Tp < Tc- In this paper we explore in some detail this phenomenon and we clarify what dynamical 
process is really described when the N oo limit is taken first. This helps to understand what correct use is to be 
made of the large- iV model in this area of non equilibrium statistical mechanics. 

The gross features of what goes on in the quenches below can be described with the help of Figure 2. The 
phase-ordering process of a system with finite A'^ is represented by path / connecting the disordered states A to 
the ordered states B. These latter states are mixtures of broken symmetry states. If, after equilibrium has been 
established, the N oo limit is taken as depicted by path //, the system is brought into a state D which is the 
mixture of the N —^ oo limits of each one of the broken symmetry states. If, instead, the N oo limit is taken 
before the quench, the process starts from a disordered state C of the system with infinitely many components and 
the ensuing dynamical process does not connect C to D along path ///. As a matter of fact, the process depicted 
by III does not exist. Rather, the dynamical evolution follows path IV leading to low temperature states E which 
are quite distinct from D. In other words, the system with N — oo supports two different low temperature phases, 
whose realization depends on the order of the limits t — > oo and N —^ oo. The distinction between these two phases is 
reminiscent of the difference between the zero field low temperature states in the spherical model [ p2| and in the mean 
spherical model jl^]. In particular, states E are very similar to the low temperature states in the ideal Bose gas, as 
it will be clarified in the next section. The point to be stressed here is that in the static 1/N expansion states A and 
B are reached, respectively, from states C and D, while 1/N corrections over states E are not informative on states 
B. This clarifies why the 1/N expansion can be used for quenches to Tc, but not below Tc, as an approximation for 
processes with finite N. 

What, then, is the use of the large- model for growth kinetics. There is an obvious intrinsic interest, once it is clear 
that although not describing a phase-ordering process of the usual type, yet the model is well defined and describes 
the relaxation across a phase transition. The growth process generated in the time evolution can be studied in detail 
and produces non trivial behavior. The outcome is quite interesting since, by modulating the initial noise and the 
final temperature, remarkable crossover phenomena are obtained. Here is where the model gives information also on 
systems with finite N, even if it is not perturbatively close to the phase-ordering processes. In fact, the phenomenology 
of the structure factor exhibits features which are also found in the preasymptotic behavior of systems with finite 
N |l^,|T^. Therefore, through the large- iV model insight can be gained into the very complex time regime preceding 
the onset of scaling in realistic systems. 

The paper is organized as follows. In Section II the model is introduced, the non commutativity of the t — + oo, 
N ~f oo limits is clarified and the nature of the low temperature phases in the large- A'^ model is investigated. In Section 
HI the numerical solution for the structure factor is presented and crossovers between different scaling behaviors are 
analyzed by means of the multiscaling analysis. Conclusions are presented in Section IV. 



II. THE LOW TEMPERATURE PHASES 



In the following we consider the relaxation dynamics of a system with an A^-component order parameter (l){x) — 
{(j)i{x), . . . , (j)N{x)) which is initially prepared in a high temperature disordered state and it is suddenly quenched to a 
lower temperature. The evolution of the order parameter is governed by the time-dependent Ginzburg-Landau model 

ot 0(p(X) 

where p = for NCOP, p ~ 2 for COP, and ry(a;, t) is the gaussian white noise with expectations 

< ■r]{x,t) > = , . 

<r]a.{x,t)r](3{x',t') > = 2Tp{iV)PSc,f3S{x-x')S{t-t'). ^^''^ 
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The free energy functional is of the form 



n[(t), N] = / (Tx 



V 



(14) 



where V is the volume of the system and r < 0, g > 0. The order parameter probability distribution in the initial 
state can be taken of the form 

Po [$, ^] = ^ exp I - ^ 1 d'^xcf^ {x) I (15) 

describing the absence of correlations at high temperature 

< (/.,(f,0)</);3(a;',0) >^ Mo.p5{x-x'). (16) 

As mentioned in the Introduction, if one wants to consider the N oo limit, in order to determine the nature of 
the final equilibrium state attention must be payed to the order in which the N ^ oo and t oo limits are taken. 
Let us consider first the sequence limAr^oo linit^oo- Keeping N finite, the equation of motion ( p^ ) induces the time 
evolution of the probability distribution from the initial form ( |l5| ) toward the Gibbs state 

P[$, t, N] — > Peg [$, ^] = I exp I - ^n{$, iV] I . (17) 

In the infinite volume limit Peq[4>,N] describes a disordered pure state if Tp is above Tc and the 0{N) symmetrical 
mixture of the broken symmetry states if Tp is below Tc. If we now take the iV — > oo limit (path of type // in Fig. |^), 
for Tp > Tc we obtain the pure phase C 



Peq[c^.^] = -exp<( >](fc2 + r + 5^)0(fc).0(-fc) ) (18) 

where S is given by the self-consistency condition 

'^=-^E<'^/3(^)</>/3(-^)> (19) 

k 

and /3 denotes the generic component. Below Tc, denoting by rfi the expectation value of 0(a;) in the broken symmetry 
state, the N oo limit D of the mixture is obtained 



Peq[0,cxD] — J drrip{rn)ii[(j) \ m] (20) 
where p(m) is the uniform probability density over the sphere of radius m and the pure state ^[(j) \ ffi] is given by 



^[(j) I m] = 5[{(j) - rfi) ■ m]— exp { > ^(fc^ + r + gm^ + gSi_)(l)±{k) ■ (l)±{-k) } (21) 

where = (j) — {(p ■ m)rh. The quantities S± and m are determined by the self-consistency relations 

S± = ^Y.< > (22) 

k 

r + g{m^ + Si_)^0. (23) 



In the end, computing averages with the weight functions ( |18D (|20D, we find the well known result of the large- iV 
model 
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for the equilibrium structure factor in states C and D. The average value of the order parameter is given by 



2 r (Tc-Tf 



with 

r{d-2) 

where A is a wave vector cutoff and = [2'^-^Tr'^/^r{d/2)]-\ The S{k) term appearing in (E4), below T^, is the 
Bragg peak due to ordering in the low temperature phase D. 

Let us now consider the limits in the opposite order. Taking the N ^ oo limit at the outset amounts to take the 



limit on the equation of motion (12) which becomes effectively linearized 

■kP[e + R{t)]^f,{k,t)+rfi,{k,t) (27) 



d(f>f){k,t) 



dt 
with 



R{t) = r + gS{t) (28) 

and 

= E < '^Z'^^' t) > . (29) 



V 

k 



Due to the linearity of Eq. (E^) , the time dependent probability distribution is gaussian 



P[0, t, ex.] = ^ exp I ^ C-\k, t)${k) ■ ${~k) I 



(30) 



with 



C{k,t)^<^p{k,t)(^p{~k,t)> . (31) 



Taking next the t oo limit 



P[lt,^] ^ Qeq[$,^] = |exp|-i^C-,i(fc)0(fc) •0(-fc)| 



(32) 



we obtain a gaussian state for any final temperature and it is legitimate to ask whether still there is a phase transition. 
From ( p7| ) it is straightforward to obtain the equation of motion for the structure factor 



dC{k,t) 
Ft 

which, after equilibration is reached, yields 



= -2kP[k^ + R{t)]C{k,t) + 2kPTp (33) 



- fcjf^ (34) 



with NCOP and 



J f C(fc = 0,i = 0), forfc-0 
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with COP, where hnit^oo Rit) = ^ is the inverse square equihbrium correlation length. From (pa) 



V ^A:2+^-2 



(36) 



where the sum extends over all values of k for NCOP and over fc ^ for COP. We now analyze Eq. ( |36| ) in the NCOP 
case referring to Appendix II for the modifications in the argument required by the conservation law. With V finite 
the solution of Eq. (^6[) for ^ is finite for any temperature. As expected, there is no phase transition in a finite volume 
system. In the infinite volume limit (pq) becomes 



= r + gTFB{r^) 



where the function B{x) defined by 



B{x) 



1 



{2TTf k'^+X 

is a monotonously decreasing function of x with a maximum value at -B(O) = K^A'^^^ / (d — 2). In writing (| 
fc = contribution to the sum (^6|) has been explicitly separated out. Defining 7^ 
temperature Tc = — which coincides with (|2^) , Eq. ( |37| ) can be rewritten as 



(37) 



(38) 



the 



and introducing the 
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with the solution (Appendix II) 



Tc-Tf 



TF[B{r^)-B{ 



(39) 



> 0, 
= 0, 

- 0, 7' = -r/g ( 



7^ = 0, 
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for Tf > Tc 
for Tf = Tc 

for Tf < Tc 



(40) 



which shows the existence of the phase transition at Tc- For the structure factor this implies 



Ceq{k) 



2>+7''5(fc) Tf<T, 



(41) 



and taking into account the form ( ^ ) of 7^, (^Tj) is identical to (|2J). Thus, as far as the structure factor is concerned, 
the same result is found irrespective of the order of the limits t — s- 00 and N 00. However, comparing the states, 
Qeq[(f>, 00] coincides with Peq[(l), 00] above Tc, but not below where 



1 



ieql 



exp 



I k=£0 



(42) 



This is a state exactly of the same form of the zero field low temperature state in the mean spherical model, while 
state D is a mixture as in the spherical model ||l^ . 

Thus, despite the formal similarity, the Bragg peaks in (^) and (^ ) have different physical meanings. In the former 
case, it signals the formation of a mixture of ordered states, while in the latter it is due to the macroscopic growth 
of the k = term in the sum (36). What we have here is a low temperature phase obtained by condensation of the 
fluctuations at /c = 0, as in the ideal Bose gas, with Ceq{k — 0) playing the role of the zero momentum occupation 
number. 

Finally, notice that with d = 2 the critical temperature (^) vanishes. Hence, all states with > are disordered 
states and the limits t ^ oo and N 00 commute for the quenches to any finite final temperature. Conversely, 
Tf = is not a critical temperature, rather it is an ordering temperature. Therefore, the quench to Tp = is an 
ordering process and the limits are not supposed to commute in this case. 
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III. SCALING BEHAVIORS 



In this section we investigate the time evohition of the structure factor by solving numerically the equation of 



motion (33). It is convenient to comment beforehand on the structure of this solution. Integrating (p3|) the structure 



factor can be written as the sum of two contributions 

C{k,t) = Ca{k,t)+CTik,t) (43) 

with 



Co(fc,t) = Aexp|-2/cP^ ds[P + i?(s)]| 



(44) 



and 

ft f ft 1 

(45) 



CT{k,t) = 2TFkP dt'exp|-2F^ ds[k^ + R{s)\^ 



The asymptotic behavior is analytically accessible and it has been derived in detail in Ref. [|10|. For quenches to 
Tp < Tc with NCOP, the large time behavior is given by 

Co{k,t) = ^L^ h{x) (46) 

CT{k,t)^TFL^Mx) (47) 



ith Lit) = i^/^ X = kL and 



I 4 - (i, for Tf = Tc (.r.^ 



Ux) = e-- (49) 
fT[x)= C dy[l-v)--'^e'^'y. (50) 



Jo 

Notice that both contributions are in the scaling form (|l|) with Cxik^t) dominating in the quenches to while 
Ca{k,t) dominates for Tp < Tc- One may go one step further regarding the order parameter as the sum of two 
contributions <j) = a -\- Q whose correlations account for the two pieces in the structure factor Cf){k,t) ~< upup >, 
Cxik, t) =< > and < apCj >— 0. Therefore, for Tp < Tc a, sort of two fluids picture of the quench is obtained, 
with the condensate a and the thermal fluctuations ^ having a distinct individuality due to the different scaling 
properties. Notice that the irrelevance of the thermal fluctuations is due to a dominating C and obeying the zero 
temperature equation of motion. 

With COP there is additional structure since the thermal contribution contains itself two different pieces 

C (k t) ^ I '^<(^'*)' ^ < x*k^{t) , 

' ' \c>(fc,t), for k > x*k,n{t) ^ ' 



with the asymptotic behaviors 



Co{k,t) = ALP'^^''^ (52) 



C<(fc,i) = (53) 



Cy{k,t)^^L^ (54) 
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for Tp = Tc 

p = { d-2, for < Ti. < Te (55) 
for Tf = 0. 

For simplicity, in writing ( p^^54| ) we have neglected the logarithmic difference between L(t) and km(t), putting x = kL 
and L{t) = t^/"^. The novelty is in the region < Tp < Tc, where there is a sharp distinction between what happens for 
X < X* and for x > x* . In the first case C^{k,t) dominates over Co{k,t), while in the second case Cy{k,t) dominates 
over Ca{k^t). Therefore, with COP we are led to regard the order parameter as the sum of three contributions 
(/i = (T + C< + C> which, again, are characterized by distinct scaling properties and whose correlations are responsible, 
respectively, for Co{k,t), C^{k,t) and Cy{k,t). The remarkable qualitative difference with the NCOP case, is that 
now the condensate also is of thermal origin because the Bragg peak is formed by C<- Hence, the temperature is 
not an irrelevant variable. Furthermore, the thermal fluctuations C> which obey standard scaling, as the time goes 
on propagate from the large wave vectors toward the small wave vectors, following a pattern which is important, as 
we shall see below, for understanding what happens in the realistic systems. Finally, the a contribution which was 
responsible for the condensate with NCOP, here is subdominant, but it can give rise to interesting preasymptotic 
behaviors if A and Tp are appropriately chosen. 

In the numerical study of the structure factor, we shall devote particular attention to the transition from preasymp- 
totic to asymptotic features. This we do both for NCOP and COP with d — 3 and d — 2. Parameters of the quench 
are the final temperature Tp and the strength A of the fluctuations in the initial state (p^. In particular, we will 
consider the two cases A = (small A) and A = —r/g (large A), where \/—r/g is the equilibrium value of the order 
parameter at zero temperature. The final temperature of the quench is important in two respects. First of all, for 
Tp > Tc the correlation length ^ is finite, while for Tp < Tc it is infinite. This is important because the general 
structure of the time evolution is determined by the relation between a microscopic length Lq and the correlation 
length ^ in the final equilibrium state. The initial fast transient, with no scaling, lasts up to some time to. At this 
point equilibrium is established over the length scale Lq and, if this is of the order of magnitude of ^, final equilibrium 
is reached over the whole system as well. Instead, if ^ ^ Lq, a second regime is entered [0, during which the 
scaling relations (1) and (2) hold. This lasts up to the time ii such that L{ti) ~ ^, when global equilibrium is again 
established. Clearly, if ^ is infinite, equilibrium is never reached and the scaling regime lasts forever. The second 
important feature involving the final temperature is that, while for Tp > Tc only the thermal fluctuations grow, for 
Tp < Tc there is growth of both the condensate and thermal fluctuations. 



A. Multiscaling analysis 

The interplay of all these elements produces a rich variety of behaviors which can be efficiently monitored through 
the multiscaling analysis. This works as follows: let us assume that the structure factor can be written in the general 
multiscaling form 

C(fc,i) = [£i(i)]"(^'F(x) (56) 

with X — kC2{t) and where Ci{t) and £2(i) are two lengths. The functions a{x) and F{x) are to be determined. In 
order to check on the assumption, the time axis is divided in intervals {ti,ti + r^) which in practice may also be of 
variable length n, and within each interval the logarithm of C{k, t) is plotted against the logarithm of Ci{t) for a fixed 
value of X. By measuring the slope and the intercept of the plot a{x, ti) and F{x,ti) are obtained. The procedure is 
then repeated for different values of x and over different time intervals. If a{x,ti) and F{x,ti) do not depend on ti 
the assumption ( ^6| ) is correct and scaling holds. Specifically, standard scaling is the case where Ci{t) — £2(0 and 
a{x) does not depend on x. In the case of multiscaling a{x) does depend on x and the two lengths Ci(t) and C2(t) 
differ by a logarithmic factor. Conversely, if a{x,ti) and F(x,ti) do depend on ti, scaling does not hold. Notice that 
in the case that equilibrium has been reached, the disappearance of the time dependence shows up as a{x, ti) = 0. In 
the following we shall not be interested in the determination of F(x) and we shall concentrate on a{x). 

As an example, consider the quench to Tp = 0. With NCOP the exact form of the structure factor is given by Q 

C(fc,t) = Aexp{-[Q(i) + (fcL)2]} (57) 
where Q{t) is a function of time and L(t) = t^^^. Hence, we have the natural choice £2(0 = L(t) and, using (p6|). 
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which shows that if there is scahng it is of the standard type. This occurs for large time where Q{t) = —d\ogL{t) 
suggests to take Ci{t) = C2{t) = L{t), yielding a{x,t) = d. Conversely, for short time, if A is small enough to allow 
for the application of the linear approximation, we have Q{t) = 2Tt and there is no scaling since the time dependence 
does not drop out 

With COP, instead, the exact form of the structure factor is given by ^ 

C{k, t) = A cxp (fc„L)V(fc/fcm) (60) 
where km{t) is the peak wave vector, L{t) — i^/* and ^ is given by (||). Choosing C2{t) — k:^{t) we find 

_ [k„.{t)L{t)]' i^jx) 

logA(t) > 

which shows that the x dependence cannot be eliminated and therefore that scaling can only be of the multiscaling 
type. This occurs for large time with Li{t) = (k'^'^L'^^^^'^ yielding 



Using (fTl), the two lengths Ci{t) and C2{t) are related by ^ (logL)^/"*. Conversely, in the very early stage 



a[x,t) = dij{x). (62) 

£2(t) 

where the linear approximation holds C2{t) = k^ ~ yJ—2/r is time independent and Ci ~ t^^'^'^ yielding 

(^{x,t) -!—ip{x) (63) 
logi 

which displays the absence of scaling through a time dependent prefactor in front of ijjix). 



B. Evolution of a{x,t) 

In the following we illustrate the evolution of oi(x^t) obtained numerically over a sequence of time intervals, with 
r = — 1 and g — 1- 

NCOP d = 3 

Let us begin with NCOP in the three dimensions. The critical temperature is finite Tc = 2tt, the exponent a is 
obtained setting Ci{t) = C2{t) and extracting this length from the inverse of the halfwidth of the structure factor. 
According to the general outline presented at the beginning of this section, we expect to detect the establishment of 
equilibrium for Tp > Tc through the vanishing of a and the scaling behavior lasting indefinitely for Tp < Tc, through 
the disappearance of the time dependence around a non vanishing value of a. This is clearly illustrated in Fig. 
where a is plotted for a quench well above Tc with Tp = 20 and for a quench to Tc- In the first case the lines rapidly 
collapse on a = 0, while in the second case the collapse is on a = 2. It is then interesting to consider the case of 
a final temperature slightly above Tc, which corresponds to a large but finite ^. Also in this case a is expected to 
collapse eventually on a = 0, however as long as L{t) is large but smaller than ^, one expects to observe a behavior 
similar to the one in the quench to Tc. Indeed, this is what happens in Fig. ^ (panel a), obtained by plotting a for 
Tp — 6.35. After the initial transient there is a collapse on a = 2, revealing critical scaling. However, this does not 
last indefinitely as in the case of the quench to Tc, but it lasts for the time necessary for L{t) to catch up with ^. 
After this a new transient sets in and the eventual collapse on the equilibrium value a = takes place. For Tp < Tc 
the behavior of a is quite similar to the one for Tp — Tc, since in both cases ^ = 00. The only difference is that the 
asymptotic behavior produces collapse on the value a — d — 3. For Tp slightly below Tc, e.g Tp = 6.20 in Fig. ^ 
(panel b), there is crossover from critical scaling with a = 2 to the final value with a = 3. In general, the asymptotic 
behavior for < Tp < Tc and Tp = is the same (Fig. ||), confirming the irrelevance of thermal fluctuations. In all 
quenches considered, the variation of the size A of initial fluctuations does not produce significant differences. 

NCOP d = 2 
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In two dimensions the critical temperature ( p6[ ) vanishes. So, for any Tp > one should observe a behavior for a 
similar to the one obtained with d = 3 and Tp > Tc- This is the case for quenches with a final temperature well above 
zero, where behaviors of a very close to the one in panel a of Fig. || are obtained. Similarly, for the quench to = 
the same behavior of Fig. |5| (panel b) is obtained, except that now the lines collapse on a = d = 2. 

A case to be considered separately is when Tp is finite but very close to zero. Fig. ^ corresponding to Tp — 0.3 and 
Tp = 10^®, displays an intermediate scaling behavior with a — 2 preceding the eventual collapse on the equilibrium 
value a = 0. At first sight this looks like the behavior of Fig. |^ (panel a) in the quench to a final temperature 
slightly above Tc. However, the interpretation is more subtle, since what is growing here, as long as L{t) < ^, are not 
the critical fluctuations but the condensate. The distinction between condensation and critical behavior is actually 
impossible to make with NCOP on the basis of the value of a, since with d = 2 in both cases a = 2. This remark will 
become clear with COP, because in that case the growth of the critical fluctuations is associated to standard scaling, 
while the growth of the condensate gives rise to multiscaling. 

COP d = 3 

We now move on to the case of COP with both small and large initial fluctuations. As discussed above, a(x, t) 
is extracted by defining x = k/km, where fcm(i) is the peak wave vector, and plotting logC vs. log£i(i) with 
Ci{t) = {k'^'^L'^Y^'^ and L{t) = t^/^. For quenches to Tp > Tc the behavior of a{x,t) essentially follows the same 
pattern as in the NCOP case. Apart from some differences in the time dependent transients, again for Tp > Tc and 
for Tp — Tc curves collapse, respectively, on a = and a = 2. In Fig. ^ we illustrate what happens in the quench 
slightly above (panel a) and slightly below (panel b) Tc- In both cases there is a preasymptotic standard scaling 
behavior with a — 2 due to the critical point in the neighborhood, followed by the incipient crossover toward the 
asymptotic behavior. Above Tc this is a = 0, as revealed by a deviating downward, while below the deviation 
occurs upwardly, for x < x* , toward the asymptotic form of Fig. 

Where the difference between NCOP and COP becomes remarkably evident is in the quenches well below Tc- Let 
us first consider (Fig. ||) Tp = 0. After a time dependent transient (panel a), which for small A in the early stage is 
well described by (|6^), the curves of a{x,t) collapse (panel b) on the limiting curve ( |62| ) depicted in Fig. |l|. Instead, 
for Tp — 1 (Fig. the collapse occurs on the finite temperature asymptotic form of Fig. ^ with minor differences in 
the transient due to the size of A. However, going to a temperature much lower but finite (Fig. |l^), while for A = 
(panel a) a follows the same pattern as in Fig. ^, the behavior in panel b with A = 1 is drastically different. What 
we have here is multiscaling as in the quench to Tp = 0, for x < x* , and standard scaling with a = 2 for x > x* . 
All these features can be accounted for on the basis of the discussion made at the beginning of this section. As long 
as Tp is sufficiently large or A small, as in Fig. ^ and in panel a of Fig. ^ only C< (fc, t) and C> (fc, t) contribute to 
a{x, t) producing the characteristic behavior of Fig. |^. However, when A is finite and Tp small enough, there can be 
a sizable interval of time during which C^{k^t) dominates over C^{k,t), for x < x* , producing the scaling pattern of 
panel b in Fig. ^ This behavior is preasymptotic and eventually the crossover to the pattern of Fig. || takes place, as 
C< (fc, t) grows large enough to overtake Co(fc, t). In our numerical solution the computation was not run long enough 
to actually detect this crossover, but it is clear that by modulating the parameters of the quench, A and Tp, the 
crossover time can be varied at will. 

COP d=2 

It is now interesting to see how this variety of behaviors is affected by pushing the critical temperature to zero in two 
dimensions. The novelty with respect to the previous case is that now the line of fixed points in between Tp = and 
Tc has disappeared and with it, supposedly, also the associated asymptotic behavior. Actually, only the fixed point 
at Tp — has survived. Thus we should observe either the relaxation to equilibrium for Tp > 0, or the multiscaling 
behavior for Tp = 0. Indeed, for temperatures like Tp = 10 the collapse on a = is observed. Similarly, the quench 
to Tp = produces a behavior identical to the one in Fig. ^, except that the peak value of a is given by 2 in place 
of 3. The interesting novelties arise when quenches to small, but finite, Tp are considered. In fact, when Tp is low 
^ is large and an intermediate scaling regime preceding the final equilibration is expected. The question is, what it 
will be like. For Tp = 10^^ (Fig. |ll|) and A = there is standard scaling with a = 2, while for A = 1 the pattern is 
identical to the one in panel b of Fig. |lO|, except for the obvious modification a{x = 1) = d = 2. The first observation 
is that these are preasymptotic behaviors, since eventually a must vanish. The second is that this phenomenology 
can be understood regarding, as long as L{t) < ^, the structure factor as made up of the three contributions (|5^^. 
Namely, as long as L{t) < ^, there is growth of the condensate and of thermal fiuctuations, as in the quench below the 
critical point. Thermal fiuctuations are clearly due to Tp > 0, while the growth of the condensate originates from the 
underlying fixed point being at = 0. The net result is that the intermediate scaling behavior with A = is due to 
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C<^{k,t) and Cy{k,t), which scale both Uke since p vanishes when d = 2. Instead, with A = 1 we are confronted 
again with a situation where, in the time of the computation, Co{k,t) dominates over C^{k,t) producing the pattern 
of panel b. As a matter of fact this is a pre-preasymptotic scaling behavior, since when C<c{k,t) overtakes Co{k,t) a 
behavior of the type in panel a is expected to occur before the eventual relaxation to a vanishing a. 

IV. CONCLUSIONS 

In this paper we have investigated the origin of the non commutativity of the limits i — > oo and ^ cxd in the 
dynamics of the first order transitions. The main result is that when the N —^ oo limit is taken first the underlying 
phase transition, which we have called condensation, is qualitatively different from the usual process of ordering 
obtained with the limits in the opposite order. In particular, condensation in conjunction with COP gives rise, 
for reasons which are not yet clear, to two phenomena which are strikingly different from what one has in phase- 
ordering, namely i) multiscaling and ii) relevance of the thermal fluctuations. We have then proceeded to an extensive 
investigation of the scaling properties of the structure factor in the large- iV model for quenches to a final temperature 
greater, equal or lower than T^. We have found a rich variety of behaviors, which can be studied in great detail through 
the multiscaling analysis. Particularly interesting is the existence of preasymptotic scaling which can be explained 
through the competition between different components of the order parameter with distinct scaling properties. 

Even though it is quite clear that the large- A'^ model is not perturbatively close to the phase-ordering processes in 
realistic systems, in concluding the paper we wish to elaborate on the connections that nonetheless exist. In order to 
do this we use the BH model [|ll| as an intermediate step. As mentioned above, in this model the structure factor for 
the quenches to Tp = with COP displays standard scaling for any finite value of A^. The discussion in the previous 
section on the behavior of a in the quenches to Tp = 10~^ helps to understand how this comes about. By integrating 
formally the equation of motion, the BH structure factor can be written as the sum of two contributions 

Cik,t) ^Coik,t) + CNik,t) (64) 

where Co{k,t) is given by while 

CNik,t) = f dt'R{t')^^^D{k,t') (65) 

Jo Co[k,t'j 

with 

D{k,t)^ J ^^C{k-h,t)C{h~k2,t)C{k2,t) (66) 

contains the nonlinearity and R{t) is given by (^8|) and (^9|). Although Eq. ( p3| ) refers to a quench to Tp > with 
N = oo and Eq. ( |64| ) to a quench to Tp = and A^ finite, the mechanism regulating the competiton between the two 
terms is the same. In particular, in both cases Co(fc,t) can compete with the second term only for x < x*. BH have 
shown that C]\[{k,t) asymptotically obeys standard scaling with a = d. It is then clear that by choosing A and A^ 
properly there may be a preasymptotic regime during which Co(fc, t) dominates for a; < a:*, much in the same way as in 
the previous section Co(fc, t) was found to dominate over Crik, t). Here 1/A^ plays a role similar to that of Tp and the 
crossover time depends on both A and A^ [ pl1Jl4| . With A^ = 100 and d — 2 (Fig. ^2|) a behavior for a{x, t) is obtained 
which is practically the same of that in Fig. |11|. In order to complete the picture, we reproduce in Fig. n3 the behavior 
of a{x^ t) for the scalar system obtained in Ref. ||l^ by the simulation of Eq. (|l2|) with A^ = 1 and d — 2. Again the 
same pattern is found, revealing that the same mechanism is operating also in this case. Therefore, one may conclude 
that the behavior of the large- A^ model associated to Co(fc, i) describes what happens in the preasymptotic regime 
also in the phase-ordering processes over the shrinking range of wave vectors with k < x*km{t). In other words, the 
asymptotic regime is preceeded by a time regime where phase-ordering over the short length scale seems to coexist 
with condensation over the large length scale. This is not surprising for NCOP since correlations are established over 
regions of size L{t), and the statistics can be expected to become gaussian over distances larger than L{t). For COP 
it is less straightforward, although the occurrence of gaussian statistics on large length scales can be detected much 
more easily through the appearence of multiscaling behavior. 
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V. APPENDIX I 



Prom the real space scaling form of the correlation function is given by 

G{T,t)=r"~''g{r/L{t)) (67) 

where g{x) is the scaling function and g(0) is a finite quantity. For r <^ L{t) ( p7| ) gives the equilibrium decay of the 
correlation function 

Ge,(f)~r"-'*. (68) 
On the other hand, the correlation function on a fractal [Q, decays as 

Ge,(r')^r2(^-'^) (69) 



where D is the fractal dimensionality. Comparing (|69|) with (J^a^ the fractal dimensionality of the correlated regions 
is given hy D — ^{a + d) and using (^) 



VI. APPENDIX II 



NCOP 

Equation (p7|) can be rewritten as 



and using r + gTcB{0) = 



where 



r'^r + gTpB{0) + + gT^ [B{C') - 3(0)] (71) 



^=c + TF[Bir^)^BiO)] (72) 
9 



Since [i?(^~^) — B{Q)\ is a non positive monotonously decreasing function, there is a positive solution of ( fz^ ) for 
c > 0, a vanishing solution for c = and no solution for c < 0. For Tp < the quantity c cannot be positive, because 
in that case would be positive and the second term in the right hand side of ( [73| ) would vanish in the infinite 
volume limit, producing c < 0. Therefore, c can only vanish implying 

,1™ TTT^ = — ^ • ( "^4) 

v^co VS, ^ g \ Tc J 

COP 

With conserved order parameter, Eq. ( ^6|) is replaced by 

which allows for a solution > —k'^in, where kmin ^ V~^l'^ is the minimum value of the wave vector. For 
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Tf < n 



-y- 

V ^ 



(76) 



the solution is negative — < ^ ^ < 0. In the infinite volume limit Tc — ^ 7c and ^ 
® as 



'Klin- Thus, rewriting 



gTp 



we can analyze this equation exactly as in the NCOP case, obtaining for Tp < Tc the analog of 



(77) 



(78) 
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FIGURE CAPTIONS 



1. Spectrum of the multiscaling exponent a{x) for quenches with COP and d = 3. 

2. Schematic representation of the relaxation processes in the systems with N finite and N = oo. 

3. Evolution of a{x) for NCOP, d = 3. (a) Tp ^ Tc (6) Tp — Tc- Different curves refer to a sequence of time 
intervals growing exponentially with the label. In this and all other figures except Fig. ^ very early times are 
not shown for simplicity. 

4. Evolution of a{x) for NCOP, d = 3 for quenches to Tp slightly above and slightly below T^- Symbols and time 
intervals are related as in Fig. ^. 

5. Evolution of a{x) for NCOP, d = 3 for quenches io Q < Tp < Tc and Tp ~ 0. Symbols and time intervals are 
related as in Fig. ||. 

6. Evolution of a{x) for NCOP, d = 2 for quenches to Tp = 0.3 and Tp = 10^®. Symbols and time intervals are 
related as in Fig. |^. 

7. Evolution of a{x) for COP, d — 3 for quenches to Tp slightly above and slightly below Tc- Symbols and time 
intervals are related as in Fig. ^. 

8. Evolution of a{x) for COP, d = 3 for a quench to Tp — 0. (a) early times (b) intermediate to late times. 
Different curves refer to a sequence of time intervals growing exponentially with the label. 

9. Evolution of a{x) for COP, d = 3, Tp = 1. (a) zero initial fluctuations, (b) large initial fluctuations. Symbols 
and time intervals are related as in Fig. ^. 

10. Evolution ofa{x) for COP, d = 3,Tp = 10"^. (a) zero initial fluctuations, (b) large initial fluctuations. Symbols 
and time intervals are related as in Fig. R. 
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11. Evolution of a{x) for COP, d = 2, Tp = 10 ^. (a) zero initial fluctuations, (6) large initial fluctuations. Symbols 
and time intervals are related as in Fig. ^. 



12. Evolution of a{x) for the solution of the Bray-Humayun model with N — 100, d — 2, Tp — 0. (a) small 
initial fluctuations, (6) large initial fluctuations. Different curves refer to a sequence of time intervals growing 
exponentially with the label. 

13. Evolution of a{x) obtained in Ref. |l4) by simulation of a system with iV = 1, d = 2 and Tp = 0. (a) small 
initial fluctuations, (6) large initial fluctuations. Different curves refer to a sequence of time intervals growing 
exponentially with the label. 



15 




N finite N infinite 



6.0 



(a) 



4.0 



T = 20 



2.0 



0-0« ^ <>^^ ^<K> <M> O-O^ 



<s^<^ o-<y^ o 



0.0 



1.0 



£r-/\-A i?r-AT\ A-A-A A A-A-& A-A-A A-A-A A 



A- 



2.0 



3.0 



6.0 



4.0 



0.0 



0.0 



T = T 



1.0 



(b) 



A- 



2.0 



3.0 



X 



6.0 



(a) 



4.0 



T = 6.35 



Eg— K3 ^■sg— Eg sj sg— sg -sg sg— sg-sg Kg-sg— sg sg- sg— fJJ -VJ^^ BiJ IKJ ^ 



- +-+-+-+ - 



0.0 



0.0 



1.0 



2.0 



3.0 



6.0 



4.0 



2.0 



0.0 



0.0 



(b) 



T = 6.20 



't+-++ 



1.0 



2.0 



3.0 



X 



6.0 



T=1 A=0 

4.0 



¥-V-V V-V-V ¥■ V-V-V V-V-V V-V-VV- V-V-V V-V-V V-V-V V- V-V-¥ V-V-V V-V-V V-V-V V-V-V '^-w-o 



2.0 



6.0 



(a) 



0.0 

0.0 1.0 2.0 3.0 



(b) 



T = A =10"' 

4.0 

'yo^ -^o^ ^^^^ ^^^^ <x>^ <y'>o <^'>o^ 'X>^ -^oo OKyo^ ^<^^ o-^^ ^^^^ 



2.0 



0.0 I 

0.0 1.0 2.0 3.0 

X 



6.0 I ' , ' , ' 

(a) 



4.0 - T = 0.3 



2.0 



0.0 

0.0 1.0 2.0 3.( 

6.0 



4.0 





1 


(b) 


T = 10"' 








? ^ ^ 





0.0 

0.0 1.0 2.0 3.1 

X 




:x) 



X 



